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EFFECTS OF SIDE-EDGE BOUNDARY CONDITIONS AND TRANSVERSE 


SHEAR STIFFNESSES ON THE FLUTTER OF ORTHOTROPIC 

PANELS IN SUPERSONIC FLOW 

By Deene J. WeicLman 
Langley Research Center 


SUMMARY 


The basic linearized equations for orthotropic panel flutter are solved by 
the Galerkin method for simply supported leading and trailing edges and various 
side-edge support conditions. In addition, allowance is made for the finite 
transverse shear stiffness of simply supported panels. These qualitative 
results show that both transverse shear stiffness and side-edge boundary con- 
ditions can be extremely important in determining the flutter boundaries for 
orthotropic panels. 


INTRODUCTION 


Panel flutter has received much attention in the past few years (refs. 1 
to 8) , the analysis of uniform isotropic rectangular simply supported panels 
receiving the most attention (refs. 1 to 4) . However, aircraft and spacecraft 
are now being constructed of various built-up panels (corrugation- stiffened 
skin panels or sandwich core panels, for example), and hence an understanding 
of the flutter behavior of such orthotropic panels is necessary. Several 
recent papers (refs. 7 and 8) have initiated investigation into orthotropic 
panel flutter. A theoretical investigation was undertaken of some effects not 
considered in these references, but thought to be of importance in the predic- 
tion of the flutter characteristics of orthotropic panels. 

As the influence of core transverse shear stiffness has been shown to be 
of importance for uniform isotropic sandwich-core panels (ref. 6), the general 
case of a panel with two arbitrary transverse shear stiffnesses and 

has been investigated. The results of the investigation reported herein indi- 1 
cate that these stiffnesses can have a large influence on the flutter phenomena 
for orthotropic panels. 

Experimental results on orthotropic panels (ref. 9) have shown that the 
panel side-edge boundary conditions (along the edges parallel to the airflow) 
have a significant effect on their flutter characteristics. In order to inves- 
tigate the importance of these boundary conditions, a Galerkin solution was 



used to obtain the flutter equation for flat rectangular orthotropic panels 
with general boundary conditions. Specific solutions are also presented for 
simply supported, clamped, and free boundary conditions along the side edges, 
and it is shown that free boundaries can significantly increase the flutter 
susceptibility of an orthotropic panel, especially where the bending stiffness 
in the crossflow direction is much higher than the bending stiffness in the 
flow direction (as is often the case). 


SYMBOLS 


^n 

a,b 

T(r) g(r) 

A qp ' B qp 
A(r),B(r) 


series coefficients appearing in expansion of w 

x- and y-dimensions of orthotropic panel (see fig. 1 ) 

flutter parameters for general orthotropic panels, defined by 
equation (5) 

flutter parameters for simply supported orthotropic panels 

including transverse shear stiffness, defined by equation (25) 


A-, (p>n), Ap(p,n),A^(p,n) coefficients in free-free mode shape expression 

(see eq. (17)) 

Bi(p ,n),Bp(p,n) coefficients in free-free mode shape expression (see 

eq. ( 17 )) 

C x transverse shearing stiffness parameter in the x-direction, 

jf2 Dx 

2 a 2 % 


D x ,Dy 




D 


*y 


F pn (y) 


^hp( x ) 


transverse shearing stiffness parameter in the y-direction. 


2 D. 


JC 


2 b 2 D Qy 


bending stiffnesses of orthotropic panel 
transverse shearing stiffnesses of orthotropic panel 
twisting stiffness of orthotropic panel 

assumed pth mode shape in the y-direction (with p half-waves) 
assumed nth mode shape in the x-direction (with n half-waves) 
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T (r)/ \ 

X qp ( m ' n ) 


cross-stream mode shape integral expression., 
'h 


fe F Pn (y) iy 


o ^ iy 

i, j,m,n,p,q,r, s mode shape numbers (number of half-waves) 


J^(cup) 


streamwise mode shape integral expression. 


a - 


r_1 [ 9mq(x) ~ Gnp(x) dx 
J 0 dx 


k 

KpC^^Cj) 

L,K 

m 

Mx.My 

M 

W x ,N y 

W x y 

q. 

R x ,Ry 

v y 

w 

x,y 

^n 


frequency factor, 

factors defined by equation ( 25 ) 

constants defined by equation (l9b) 
mass of orthotropic panel per unit area 
applied bending moments per unit length 

Mach number of airflow 

applied in-plane forces per unit length for orthotropic panel 
(positive in compression) 

in-plane shearing force per unit length 



dynamic pressure 


in-plane loading factors. 


R x a 2 Ky-a 2 

— - — and — — , respectively 

knD x 


vertical shearing force per unit length 


lateral deflection of orthotropic panel 
coordinate axes (see fig. 1) 

set of transverse shearing stiffness constants (see eq. (20)) 
nondimensional set of constants (see eq. ( 23 )) 
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p = ^/m 2 


1 


D. 


■boundary condition parameter, p. x + — 


xy 


qp 


- 2^ 

9 = ^ 


% 

A 

^x'^y 
9 = 1 + 


D. 


EL- 

^x°y 


Krone cker delta 


constants defined in equation ( 19 ) 


dynamic-pressure parameter. 


2qa^ 

pD. 


x 


critical dynamic-pressure parameter for panels with infinite 
shear stiffness 

Poisson's ratios 


to 


circular frequency 


THE INFLUENCE OF SIDE-EDGE BOUNDARY CONDITIONS 


General Galerkin Solutions for Panels Simply Supported 
at Their Leading and Trailing Edges 

The orthotropic panels under consideration have a width b (in the 
y-direction) and a length a (in the x-direction) , and are subjected to a 
supersonic airflow of Mach number M flowing over the upper surface in the 
x-direction. (See fig. 1.) These panels were assumed to be simply supported 
at both the leading and trailing edges and subjected to in-plane loadings N x 
and Ny considered positive in compression. If simple static strip-theory 
(Ackeret) aerodynamics are used and simple harmonic motion is assumed, as dis- 
cussed in reference 1, the differential equation for the lateral deflection w 
of the orthotropic panel becomes 


d^w 2 / P xy \ bK ^y 3^w dfw dfw _A_ 3w 

3? V>x ^3x 2 3 Y 2 D x3^ a 2 3X 2 a 2 Sy 2 a 3 ^ " 
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Figure 1.- Loadings and dimensions for a typical 
orthotropic panel. 


If the product li x |iy is not very small with respect to 1, this equation is 
still valid if D x and Dy are redefined by dividing them both by 1 - p x Py. 
Generally, the p.y-factor in the second term of equation (l) is neglected. 

In the appendix the solution of equation (l) is determined by the Galerkin 
method for general boundary conditions. To illustrate the influence of side- 
edge boundary conditions, however, several explicit examples are computed in 
this section. The Galerkin method is used here for generality, even though in 
a few cases an exact solution of equation (l) is possible. 

For panels simply supported at their leading and trailing edges (the 
boundaries parallel to the y-axis) the deflection must satisfy 




w] „ = 0 

_l x=0,a 

^x-0. 


= -D-y 


S^w 

chc a 




3y 2 / 


= o 


x=0,a 


( 2 ) 


To satisfy these boundary conditions, the deflection was selected as 

w = Y, ®pn Sin F pn (y) 


n 
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Only a one-term expansion In the y-direction has been assumed in this case and, 
thus, a summation with respect to p does not appear. The function *pn(y) 
has the subscript n because of possible coupling of the x and y functions 
through the particular side-edge boundary conditions used in some cases. The 
Galerkin procedure is applied by substituting equation ( 3 ) into equation (l), 

multiplying the resulting equation by sin SH2L F pm (y), and integrating over 
the area of the plate. The following equations then result: 


( 1)- r(m) m 2 — (m)^ 2A V / mn \f7 ,^+n -T^PP 

(“ - 4p - K- ♦ 7 . 2, Ssnnj-p 5 - 

/ n / x PP ' ' 


= 0 


(nj/m) 




where 


-r(r)/ . .r-1 

I (m,n) = b 
IP 


0 


f (y) ~ f (y) dy 


qm 


dy 1 " 


pn' 


jM _ - WV . u V 


qp 


- ( r) 

B qp 


-2 - 1 a 2 

^ ^2 Z(0J 


1 D y a 4 I S )(r ' r) 


t 2 I^ j (r,r) ^ D x h 4 I^ } (r,r) 


qp 


J 


(5) 


For any specific boundary conditions in the crossflow direction, then, the 
functions F pn (y) must be chosen to satisfy all the y boundary conditions 

and the integrals I^^(m,n) must be evaluated. Substituting these expressions 

into equations (4) yields a set of homogeneous linear algebraic equations for 
the coefficients ap n and the condition for a nontrivial solution is that the 

deter mi nant of the coefficients be zero. This determinant yields an algebraic 
equation that must be solved for the dynamic pressure parameter A. 

These general equations (eqs. (^)) are examined for the special case of 
only two terms in the series expansion for w. Although this two-term solution 
is not exact, it does yield the qualitative results necessary for the evalua- 
tion of the effect of crossflow boundary conditions. If only a two-term 
expansion (in which m + n is odd) is used, the condition for a nontrivial 
solution yields for the flutter parameter A 
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A = 




2 2 

n - m 


| ^ - SM)(-n* + + iiW) JE 


(°)(m, m) I^ (n,n) 


I^ 0 )(m,n) l(°)(n,m) 
PP PP 


(6) 


The critical value of A is determined hy solving for the lowest value of A 


at which the frequencies for the two modes coalesce. Thus, 


dA = 

Scjo 


since 


the frequency appears only in k , this relation is equivalent to 
yields the critical value 


= 0 | si 

a§-o) 

bk 2 ) 



” 2 4 p’ 






( 7 ) 


If the boundary conditions in the y-direction are such that the function 
Fp n (y ) is independent of the mode number n in the x -direction, the critical 

flutter parameter becomes the simple expression 





( 8 ) 
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Figure 2.- The value of m for minimum A cr as a function of A for panels simply supported 
at the leading and trailing edges (Galerkin solution using m and m + 1 half-waves). 
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If, in addition, two consecutive values of m and n are chosen (n = m + l), 
the minimum A cr is found at one of the two integer values m nearest the 

value i|-l + \j 1 + l/i"- 2Ap p j. A plot of m against Ap p is shown in fig- 
ure 2, and values of A cr for several values of m and Ap p are given in 
table I. Table I indicates that m = 1 does not give the minimum for this 
process for large negative values of Ap p . However, a difference of only 
10 percent is incurred in the range investigated by using m = 1. 

TABLE I.- VALUES OF \, r FOR A PANEL SIMPLY SUPPORTED AT THE LEADING 
AND TRAILING EDGES AND EITHER SIMPLY SUPPORTED OR CLAMPED ON 
THE SIDE EDGES (TWO^MODE GALERKIN SOLUTION) 


m,n 

Values of \. r for 

values of A pp of - 

-10 

-100 

-1000 

-10 000 

1,2 

3,7 

4.5 

5.6 

6.7 

0.822 x io3 
1.167 
1.770 
2.515 
3.487 

5.753 X 10 5 

5-733 

6.215 

6.953 

7.907 

5.507 X 10^ 

5-139 

5.096 

5.134 

5.211 

5.782 x 10 5 
5.080 
7.987 
7.952 
7.971 
7.971 





Since end fixity is to be investigated, only the modes m = 1 and n = 2 
are used (even though it is possible that other modes might be critical in cer- 
tain situations). The cases of simply supported, clamped, and free boundary 
conditions along the edges parallel to the- flow are discussed. 

Panels with the side edges simply suppo rted.- The boundary conditions in 
the y-direction are satisfied by F (y) = sin E2Z. Thus, the cross-stream 
integral expressions are independent of n and can be written as 

T (r) / , \r /2 r r /l\ 

!qp = (- 1 ) P * 

since p = q. and r is even. The parameter Ap p (for the single pth term in 
the y-direction) becomes 

*pp = ^ + p2 (9) 
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The value of >^ r is obtained by substitution of App into equation (8) and 

selection of m = 1 and n = 2. The results are the same as those of refer- 
ence 1, where the lowest term in the y-direction (p = l) was -used for a uniform 
plate. These results are shown in figure 3* where the critical flutter param- 
eter ~h cr is plotted against A^j_. 



*11 


Figure 3.- Flutter results from a two-mode Galerkin solution (lowest two 
terms) for orthotropic panels simply supported at their leading and 
trailing edges. 


Panels with side edges clamped .- In the case of clamped panels, several 
functions for Fp n (y) could be used to satisfy the clamped boundary conditions 





( 10 ) 
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The Iguchi function (ref. 10) is assumed here and is defined as 


3pn<y> - $ - 1 f ♦ - l) - 


sin HZ 

D 

Pit 


(11) 


Therefore the integrals become (see ref. ll) when q = p 


With these functions, 
considered 



( 2 ) 

qp 

( 2 ) 

<3P 


( 4 ) 

qp 

■( 4 ) 



-12 + 


n 2„2 

q rt 


1 

2 


1 

2 


if the lowest term in the y-direction with 


(q odd) 


(q even) 


(q odd) 


( q even ) 


(q odd) 


(q even) 


p = 1 is 


A 11 = - 2 -472 + Hyj (12) 

The critical values of the flutter parameter are given by equation (8) by rede- 
fining A 11 according to equation (12). The flutter boundary which has been 
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plotted on figure 3 for simply supported side edges is identical to the flutter 
boundary for these clamped side edges. 

Panels with the side edges free .- Panels having free edges in the cross- 
flow direction must have mode shapes satisfying the following boundary 
conditions: 


* - §) 


0 


y=0,b 


v iUb _D y 


c)^w I ~xy' 

^5 + l^ + 




w 


D y /dx 2 3 y 


= 0 


y=o,b 


a2 w n 2^2 

and since — — ^ = w, the function of y must satisfy the boundary 

ox^ a . 

conditions : 


where 


<3- Fpn(y) n2jt2 


dy 2 


M-xFpn(y) 


= 0 


y=0,b 


<i 3 F (y) n 2 it 2 aPpnfri 


dy- 


pn2 

"dy" 


= 0 


y=0,b 


D. 


7 = Px + 


D,: 




(13) 


In order to satisfy boundary conditions (13), the general function of 
y (*pn(y)) mus l °f necessity be a function of n (the number of waves in the 
x -direction) . The function selected herein is 


F pn(y) = cos TT + 


[i + (-D p ] [ Al + Aa(g - I) 2 + a 3 (£ - | j 


(14) 
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I 


■where 


A x (p,n) 


-p 2 a 2 1 3 a 2 A 3 (p,n) / 2 8 a 2 \ 

2 b 2 p x n 2 2 jf 2 b 2 n 2 (j. x y 3 7 n 2 it 2 b 2 7 48 a 2 y 


Ag(p,n) 


12 - 


yb 2 n 2 rt 2 '' 


2a g 
yb 2 n 2 jc 2 
a 2 


|A 5 (p,n) 


BpCp^n) 


6B 2 (p,n)a 2 { ^ 

TtfW 5 ^ 


2bV£ \ 

8a 2 / 


B 2 (p,n) 



2 2 b 2 n 2 ir 2 \ 

-7 

i P * + ^ a 2 j 

6 Px| 

f JL - 1 + 7P 2 ??* 2 ) 
yPx 12 a 2 J 


for any arbitrary A^(p,n) . For even values of p with A^(p,n) = 0, the 

slope at the free ends would, become zero, an unnatural condition; however, a 
slight curvature of the panels at the side edges is caused by using a value for 
A^(p,n) in such cases. The integrals then become 


I ( 0 ) (m,n) = i 

qp ’ 2 


qp 


2 2 
q jt^ 


B-^p^n) - q 2 « 2 - 6jB 2 (p,n) - 


2 2 
p^rt 


B-j^q^m) 


Bi^m) Bg(p,n) Bg(q,m) B x (p,n) k (3 2 2 

20 + 20 P n 



B 2 (q,m) B 2 (p,n) i . , 

+ jj2 + j B 1 (q,m) Bp^n) 


(q and p both odd) 
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(o)/ v 1 4 /q 2 * 2 \ / X 4 /p 2 jt 2 \ / x ^ A 2 (-P' n ) 

Iqp'f”.-) * 2 5 qp + + ' 12 j A > (<1 ' ,n) + Vg~ 


A 5 (p,n) + A^(q_,m) 4 A 2 (q,m) A 5 (p,n) A 3 (q,m) + A^q^m) A 5 (p,n) 

57^ + p2 n 2 + 


112 


Ag(cL,m) Ag(p,n) 

+ go — + ^ A i(^ m ) A i(p> n ) 


A^Q-jm) A^(p,n) + A^Pjn) A^q^m) 
20 

A 2 (p,n) A^q^m) + A 1 (p,n) A 2 (q,m) 


(q and p both even) 


(0). A 3 (0,n) + A 5 (0,m) A 3 (0,n) A 3 (0,m) 

Iqp ( m > n ) = 1 + 55“^ + 5^ + 2^(0^) + A-L(0, m )] 

AgCOjin) + A-^O^n) A^O^n) A 3 (0,m) + A 1 (0,m) A 3 (0,n) 

+ '6 “ + ‘ 20 ~ 

A 5 (0,m) AgCo^n) + A 3 (0,n) A 2 (0,m) 

+ — — + 4 A 1 ( 0 ,m) A^O^n) 

AgCo^m) A 2 (0,n) A^Co^m) A 2 (0,n) + A^O^n) A 2 (0,m) 


(q and p both zero) 


I qp^( m > n ) = 0 


(For other combinations of q and p) 
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t / \ \ / \ ^f-Bp(pjn) 24Bp(c j,m) 

r 9P (m ' n) " ' V 5 <JP + - - l a - — g g’ ' ♦ JE 2 ( 5 ,m) 

+ 2B 1 (q,m) Bg(p,n) + Bg(q,m) Bg(p,n) (q and p both odd) 


t ( 2), , p 2 rt 2 

I aT . (m,n) = - - 8 n _, + 48 


qp 


2 qp 


A^Cq^m) Aj(p,n) 


+ 4A 1 (q J m) A 5 (p,n) 


p 2 ^ 2 q 2 rr 2 

+ 10 A^(q,m) + jj^ Aj(p,n) - ^A^q^m) 

+ [f A 3^' n ) ” + ^(p^^q,*) + JL A2(q, m )J 

(q and p both even) 

( 2 ) f~ 

IqP (m ' n ) = [l + 2A 1 (0,m) + ± ^(0,®) + ^ A 3 (0,m)| 

+ 4Ag(0,n)jl + — A 3 (0,m) + 2A 1 (o,m) + 

(q and p both zero) 

t ( 2) / , _ _ 

qP = 0 (For other combinations of q and p) 


(4) 4 4 

Iqp (m,n) = 2JL 8qp - 4p 2 rt 2 B 1 (q,m) - (3p 2 n 2 - 24)B 2 (q,m) 


(q and p both odd) 


(4) 4 4 r _ 

I qp = 8 qp + 2A 5 (q,m) p 2 * 2 - 24 + | A 5 (p,n)l 


4A 2 (q,m)^ 2 Jt 2 + 2A^(p,n)j + 96A ] _(q,ni) A (p,n) 


(q and p both even) 
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JM, I 

I qp ( m ; n ) 


48 + | A 5 (0,m) A 5 (0,n) + 96A 5 (0,n) A^Ojm) + 8A 2 (0,m) A 5 (0,n) 

(q and p both zero) 


T 00, % 

I qp ( m > n ) 


= 0 


(For other combinations of q and p) 


where 


'qp 0 

(p t q) 

’qp = 1 

(p = q) 


These integrals are seen to be functions of the mode numbers m,n of the two 
terms in the x-direction. In all factors except the aerodynamic loading 
factor, m equals n; whereas in the aerodynamic loading factor m is not 
equal to n, since an odd derivative appears only in that term. Thus, the dis- 


tinction between m and n need only be made for the expression 


I^°^(m,n) 
qp v ' 


that multiplies the A-faetor. If a two-mode solution is considered, the 
dynamic -pres sure parameter can be written, in general, as given in equation (6) 
(if m + n is odd), and the critical value of the flutter parameter becomes 
either 


A, 


cr 



(m 2 - n 2 )| 

f k 

- n 

> ♦ ^ # ♦ $?) 

8mn , 

\ 

M p^( m , n ) 

1 4° J ( n,n) 


or upon expanding 

A, 


?t^(m 2 - n 2 ) 


'cr 


|li.H^(m,n) li2\n,m) 


8 i ^ 


) 


(m > n) 


(15) 


■A - n 4 + n 2 )^ 


- **£ - 



Ipp^(n, n ) 

n 2 b 2 



■( 2 ) 


-i 


(m,m) 


pp v * ’ 


a 


"jr^b^ D x 



4v } m 

D x 

lt°)(n,n) 
[_PP x ' ’ 


PP 



(m > n) 


(16) 
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Rewriting the constants A-^, k^, B-^ , and B 2 appearing in the integral 

expressions in terms of only two parameters yields 


where 


and 


Ai(p 


\ p 2 n 2 1 , 3A 3 (p,n)^ 2 - 8 n 2 0cp 

>*) = ~ *~2= ~ p + 2^1? 9 ' 1 "PS ~hQ~ 

2n 0 2 n 0 cp \ 5 n20 


AgCp^n) = f-pr - 

\n^ 0 cp 

- B 2 (p,n)/- i + 

\ ^ n^ 0 Cp y 


B x (p,n) 

B 2 (p,n) 


cp/ p 2 jt 2 + n 2 0 

n 2 0 (p 
\cp - 1 + -j -2 


0 = 


A 2 

^ a 2 


Bv 


9 = 2 - - ! + 




} (17) 


If the zeroeth term in the y-direction (p = 0) and the lowest two terms in the 
x-direction are considered, A cr is given by either equation ( 15 ) or (l 6 ). 

As an example of the effect of free edges on the flutter boundary, the 
flutter conditions for the following set of panels (typical of corrugation- 
stiffened panels with stiffeners running perpendicular to the airflow direc- 
tion) will be computed: 

R x = N y = 0 

0 = 0.01 


cp = 10 000 


A 5 (p,n) = 1 
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II 


then 


where 


\ r = 18.264- 


15 - ( 7.147 X 10 _5 )a ( 


T (°)/ , 
.,\ 7 (r) I 00 (r ' r) 


00 


I^(r,r) 

00 


— (r) 1 00 >(r ' r) 

400 1 ( 2)7 : 

T 00 (r ’ r) 


2 a£/^y 

^2 h2\^r + 



(18) 


A comparison of the critical dynamic -pres sure parameter for panels with 

free side edges (eq. (l8)) with the corresponding expressions for simply sup- 
ported side edges (using eq. (9)) and. clamped side edges (using eq. (12)) is 
shown in figure 4. The lowest two terms in the x-direction and the lowest 



Figure 4.- The influence of side-edge boundary conditions on the flutter of 
orthotropic panels simply supported at their leading and trailing edges 
(two-term Galerkin solution). N x = N y =0. First two terms in the 
x-direction; lowest pth term in the y-direction. 


17 


I 


allowable term in the y-direction have been used in all cases. The allowance 
of a rigid-body or "zeroeth" mode shape for the free edges is the governing 
condition and can cause the lowering of the flutter dynamic pressure q by as 
much as an order of magnitude for the example series shown in figure 4. Some 

of this difference might be caused by the low value of ^or the relatively 
large value of selected for this example, but free side edges can cause 

significant flutter differences as shown. 


Four -Mode Galerkin Solutions for Panels Simply Supported 

in the Flow Direction 

The general Galerkin equations (eqs. (4)) can be examined not only for the 
two-term expansion but for the special case of a four-term expansion as well. 
Even though this solution is not exact, it yields more accurate results than 
the two-term expansion shown in the previous section. 

For simply supported and clamped side-edge boundary conditions, a compari- 
son between the two- and four -mode solutions and an "exact" solution (presented 
later) indicates that convergence is still not accomplished, and some caution 
should be exercised in using these values for design. For the free side-edge 
boundary conditions, however, the results of this analysis may be more useful. 

By still selecting only a single p term in the y-direction, and four 
terms (m, m + 1, n, n + l) in the x-direction, the flutter parameter becomes 




K -1 + 
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T W r tn T 4n+l T ln+l j 
LK 2 / 


where 


( 19 a) 


I^(n,n+l) I^°^(n+l,n) 
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I pp )(m ’ n) I pp )(n,m) I pp )(m+1,n+1) I pp )(n+1,m+;L) 


(n 2 . m s ) 2 \(n + if - (m + l)^ I pp )(m ’ m) I pp ) ( n ' n > I p£ ) ( m+1 > m + 1 ) I pp )(n+1 ' n+l) 


( 0 ), 


r(0)f 


Ipp^Kn) I^(n+l,m+l) I^(n,n+1) + l£°\m+l,n+l) l^^n+bn) I^°^(n,m) I^(m,m+1) 


(q 2 - m 2 )[(m + l ) 2 - ufj [(n + l ) 2 - n^] {(n + if - (m + l)f]l^(m+l,m+l) I^(n+l,n+l) 


PP 


PP 


PP 

(m + n odd) 


(Equation continued on next page) 
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L = 


4p^ n+1 >“) ; 4p^ m ^ n+1 ) 4p^ n,m+1 ^ I pp^ irL+1 -’ Q ^ 


jjn + l) 2 - m 2 J 2 [n 2 - (m + 1)0 2 I^(n+l,n+l) I^°^(n,n) l£°\m+l,m+l) 


PP 


Ipp^(in,m+1) I^°^(m+l,m) l£°^(n,n+l) I^°^(n+l,n) 


[]m + l) 2 - m 2 ] (jn + l) 2 - nf] (m,m) ^p (m+l,m+l) 4 p (n,n) l£ p (n+bn+l) 

Ipp^m+l^m) I^°^(m,n+1) I^°^(u,m+l) I^(n+l,n) + l£°\m,m+l) l£°\n+l,m) I^°^(m+l,n) I^°^(n,n+l) 


[(m + l) 2 - m 2 3 [(n + l) 2 - n0 (n 2 - (m + l) 2 0 [(n + l) 2 - m0 Ip p (m,m) Ip p J (m+l,m+l) Ip p (n,n) Ip p (n+l,n+l) 

(m + n even) 


— _ 1 ) ^m^m+l 


I^°^(n,n+1) I^°\n+l,n) 
PP v ; ' PP ' 9 


T lm T ln 


I^°^(n+l,m+l) I^^(m+l,n+l) 
PP ' ' ' PP v * 


2L ]j]n + l) 2 - n 2 J 2 Ipp\ n ^ n ) Ipp^ (n+l^n+l) [( n + l) 2 - (m + l)f ] 2 I^\n+l,n+l) (m+l,m+l) 


^n^n+l 


I^^(m+l,m) I^\m,m+l) ^ lT1 I^ 0 \m,n) I^°^(n,m) I 

PP PP T ‘m+l T *n+l pp ' * ' PP ^ ' 


[(m + l) 2 - m^j 2 I^(m+l,m+l) I^(m,m) (n 2 - m 2 ) 2 I^(m,m) I p °^(n,n)| 


(m + n odd) 


k = 4< 


^WWfl 


I^^(n,n+l) (n+l,n) 

PP PP V *_ 


"V^n+1 


I^°\m+l,n) I^^(n,m+l) 
PP v * ' PP v ' 


[(n + l) 2 - n 2 ] 2 Ipp^(n,n) I^(n+l,n+l) [(m + l) 2 - n 2 ] 2 l£°\m+l,m+l) l£°^(n,n) 


T) *1 
•n 'n 


PP 

I^°^(m+l,m) I^°^(m,m+l) 
PP v PP v 


r ‘m+l T 'n 


PP 

( 0 ), 


I'“ ; (n+l,m) I^°^(m,n+1) 
PP PP 


[(m + l) 2 - m 2 ]] 2 I^°^(m+l,m+l) (jn + l) 2 - m 2 ]] 2 I^°^(n+l,n+l) 

*“ -* PP PP U ~ -J PP PP 


(m + n even) 


^ = ±2 - 4v 



If, for some particular boundary conditions under investigation, the 
y-expression Fp n (y) is not an explicit function of the x-direction 

number n, the integral expressions I' ' (m, n) are also independent 


and n, and all the integral ratios used in equations ( 19 ) become 1. 
example, the expression for L then becomes simply 


mode 
of m 
For 


>(191) 
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(m + n even) 


L = 


{ 


(2m + X)(2n + 1) [ n 2 . (m + l} 2] Q n + 1)2 


7 ^} 


L = 


(2m + l)(2n + l) jj' n + j_)2 _ ( m + l)^J(n^ - m^J 


(m + n odd) 


and A 


(r) 


PP 


and 

PP 


become Ap p and Bp p , respectively. Since Fp n (y) 


is 


often independent of n (for example, simply supported or clamped side-edge 
boundary conditions), this case is discussed in some detail below. 



For a given value of A the condition = 0 determines the value 

_ cffipp _ 

of 'A cr - A plot of A against Bp p for a selected value of Ap p must be 



made, from which the lowest value of Bp p can be found at which ^— = 0 . This 
value of Bp p is substituted back into equations ( 19 ) to determine for 

that value of Ap p . If the problem is to be calculated on a computer, it is 

f m 2 ^ - 

helpful to start Bpp at the value 


n^ (nn + or) „ ^ ^ 

2 - A u App and decrease Bpp 


by Ap p X Kj until the value of A starts to decrease. All that needs to 
be known are the assumed mode shapes for particular cases. 


The lowest two terms in the x-direction for a two-mode solution were 
selected, _and the resulting critical flutter parameter was plotted as a func- 
tion of App in figure 5* Also the lowest four terms in the x-direction were 

used in a four -mode solution as noted in figure 5* The "exact" solution curve 
is the value to which the Galerkin solution would limit if the number of 
x-terms became very large. This solution is determined from the solution of 
reference 1 for simply supported plates and is based on the fact that the 
flutter boundary for panels with cross-stream boundary conditions other than 
simply supported can be computed from simply supported results by merely rede- 
fining App and Bp p . Thus, these exact results are not strictly exact, since 

only a single approximate mode shape term was used in the y-direction. It can 
be seen from figure 5 that the two- and four-term solutions agree fairly well 
with each other, but a large discrepancy exists between these solutions and the 
exact solution and indicates that these solutions are not well converged with 
only four terms. 
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THE IHFLUEKCE OF TRANSVERSE SHEARING STIFFNESSES 


Basic Equations - General Case of Simply Supported Panels 

The two transverse shearing stiffnesses Dq and Dq^ of the panels 

shown in figure 1 can he included in the differential equation for deflection, 
provided the equation is expanded to the sixth-order (see ref. 12) equation: 


0,1 5x 6 + ^ ax 4 dy 2 + ^ chc 2 ^ 


rr + a 4 " D x rij; ~ 2 j^xyl 1 " ^•v) + 2 MxI V 


^ - <xc/n„ + n _il 


N, r ■%*-- + 2N-, 


+ i D _u_w n _ qu- N„ + w + 2JU "" 

2 ^ Jb x 2 ^y 2 ^ Sy 4 Sx 6 y Sx^dy 2 ^ bx^ 


- N X 


Sx^Sy 2 


H. + K„ h * 2»„ 




+ a« N. 


♦ Ny + V =9 k ^~s + *V + “• 


chc 2 ^ 2 






-(l-^)(N x g + N y ^ + 


dx 3y 


2<1 1" c*5 w 


+ rr^ <^W , CpV 

% SxJay 2 ^ cte Sy* 1 


-aa-r-aj 


dx ^y 2 


(i - ^My) ^ 


2!” o 4 v 


+ m “ “5 S’? + 06 + & 




where N and N„ are in-plane loads that are positive in compression, and 
x y 

are functions of the shear stiffnesses; specifically. 
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ai = 2 


1 fixy^x 


% 


D % . “10 


a 3 = a l 


St + fiQ. 


_ 1 D xyDy 
4 2 DQ X 


„ _ a l 

“5 - Do 


“6 = 


a io 

^x% 


a.^ 




Dl 


d Qx 


Dv 


“9 ' I ig^ 1 - ^) + 


a 10 = DxSy " | D xy( D xMy + DyMx) 


The shear stiffnesses can he determined by a method of approach similar to the 
method shown in reference 13. The boundary conditions for simple supports (in 
which all points along an edge are restricted from moving parallel to the edge) 
are 
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wj „ = 0 


»x=0,a 


SSL. 

d Qv 


= 0 


x=0,a 


D x /cj 2 w cl 2 w 1 ^y\ 

■^Jx=o,a “ ■ rr ''y aj2 ' d Qx sr - d q w 


= 0 


x=0,a 


> ( 21 ) 


wj. ^ ,_ = 0 


]y=0,b 


Q* 


D Qxj y= o^b 

\lyO/t 


= o 


D, : 


(c£v t 5 2 w 1 5Q y ^x 


- dx 2 D Q y ^ D Q x Bx 


= 0 


|y=o,t> 


These boundary conditions are satisfied by the assumptions 

. V V — ritx . STty 

Qx = 2, L ars cos “ sin ~ 

r s 

_ V V = • nutx njty 

Qy = 2, Z, ^ sin — cos “b 

m n 

(where a rs and ^ are constants) and finally 


w = 


Application of the Galerkin 
equations 


II 


i 3 


, utx . jTty 

a-t s sin sm — L - 

a b 


(22) 


method yields the following set of homogeneous 


2k 




8q. '\ -1 i ra is , i 2 s 2 jt^ , s^jA i 2 :t 2 s^-nf 

* 4. — + + + ^ — + ^ — 

i+r odd 


+ (l - MxMy) ~ N xy~ ^ X 

i ,i 


ljrsa-Lj 


(i 2 - r 2 )(j 2 - s 2 y 


i+r odd j+s odd 


, j it , i 2 jt 2 j 2 :jt 2 , , \ 

“r ~nr “e -~ 2 ~ * “9 -~ 2 ~ 1 - - 0 


where r and s take on all values. If the usual case where =0 is 

considered and the constants are redefined as 
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a^tr 2 

OQtr 2 

“1 

d 

tl 

“8" a2 


agir 2 

_ Oo* 2 

a 2 

Ac 

CM 

rQ 

II 


a^it 2 

- H x a2 

% 


^ “x* 2 


a,, tt 2 

If a 2 

\ 

% 2 4< 

V ^ 

% 

jA 

= a 5 ~ 

v ar 

A = 2<laJ 

P\ 


a6^ 

T2 mo^a^ 

°6 

a 2 ! 2 

k — IT 


it 



S- 

II 





the following set of homogeneous equations appears: 


i 6- 4 2- 2 4 a 2 - 6 a^ - 4 „_2 2 a 2 

a-rs S- r a l - r 3 “2 - r 3 — a 3 - s — a 4 - r - 2 ^ s — 


(l - I^My) 


V + IV + 


- 

^ D x 


%(r 6R x + r4s2 p Sy) + %( r4s2R x + ^ j| Sy) + ^7 (^ 3 \c + s6 j *sr) + “sfV + ^ j R y, 


a^A 2 ^ + S 1 *' Sg- Eyj + (l - Hjjtiy^r 2 ^ + s 2 ^ Ryj + k 2 ^^ + r 2 s 2 5g + + i^ag + s 2 &p 


+ 1 - 


hcM 


4k V ira is 7-4-22-4-2-2 \ 

-j l — — g ( “ 5 1 + “e 1 s + “ 7 s + “e 1 + “96 + 1 - MxMyj = 0 


(24) 


1 

i+r odd 
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where r and s take on all values . To determine the effect of the shear 
stiffnesses , a two-mode solution was completed and compared with the two-mode 
solution with infinite shear stiffnesses. Selecting a single term s in the 
y-direction and any two terms i,r in the x-direction (one even and one odd), 
and setting the determinant of equations (2^) equal to zero yields the fol- 
lowing expression for the flutter parameter A: 


| ±2 r 2| j l ~ K 2< r ) + r2A ( r ) + B(r) y K 2 (i) - i 2 A(i) - B(i) ^ ^ 


r 2 [I + K x (rg 


i 2 |l + K x (i)] 


where 


K n 


(j) = + J 2 s 2 ag + s^ay + j 2 a 2 + s 2 Oq - 


a^ + sc^tjag+sc^ - H x |iy 


v r . 6 - , A 2 - , .2 4 ar - o a - . . 

K 2 U J^ja-L+jsog + js — a^ + s + j 


k 


Kn . [i ♦ k i(j ^ . 2s 2 £[^ (1 - ^ 


2Dx 


B(J) = [H'K 1 CjJ|(s 2 ^I V + k 2 )- 


The critical value of A is found by setting the partial derivative of A 

O p 

with respect to of equal to zero and solving for the value of of at A c 
and then resubstituting this value into equation ( 25 ) to yield 


. rth-2 pi Tr 2 A(r) - K 2 (r) K 2 (i) - i 2 A(i) 

cr “ 8ir I "" r J 1 + K-^r) + 1 + K-^i) 

V. 


Ki(r) - Kx(i) Dy s k a k 


|l + K 1 ( i )j [l + K^r)] Dx: b^ 


(i + r odd) ( 26 ) 
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The value of 7^. r for the case of infinite shear stiffnesses is 




*Hi g _ r 2)%2 + r S _t) 

8ir 


(i + r odd) (27) 


where A(i) and A(r) both become 


A - R* “ 23 





and the assumption I^My « 1 and the fact |j.yD x = p x Dy k ave Doth been used. 

To illustrate the influence these transverse shear stiffnesses can have 
on panel flutter, the flutter conditions for the following series of panels are 
computed: 


N x = 0 



= 1,000 

Dx 

|i y = 0.15 

i = 1 
r = 2 
s = 1 


These conditions correspond to a square panel (typical of the corrugation- 

^cr 

stiffened panels of ref. 7) without in-plane loading. A plot of 


* 2 D x 


against A for several values of C v f C„ = — - — 1 and C 

\ 2a 2 D o / 



is shown by the solid curves in figure 6, and the values used in this plot are 
also presented in table II. 
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(a) Infinite Dq^ stiffness. Cy = 0. 

Figure 6.- The influence of transverse shear stiffnesses on the flutter of repre- 
sentative orthotropic panels simply supported on all sides. Two-mode Galerkin 


solution: 


lOOOj i,r,s = 1,2,1 


= 0.15. 


From figure 6 several findings are evident. For very small shearing 
stiffnesses (large C x values), the flutter parameter can drop to less than 

1 percent of the value for infinite shear stiffnesses and, thus, can have an 
extremely important role in panel flutter. It can also he seen hy comparing 
figures 6(a) and 6(b) that, for large values of C x , the influence of Cy is 
negligible, particularly for small negative values of A. Thus, C x is con- 
sidered to be of more importance than Cy, and a simplified analysis including 
only C x is presented in the next section of this report. 
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(b) Small Dq^ stiffness. Cy = 1 x lO - '' - . 
Figure 6.- Concluded. 


Infinite Shear Stiffness, Dq^ 

It can be seen from figure 6 that C x is generally of more importance 
than Cy (at least for nearly square panels of this stiffness ratio). By con- 
sidering the Dq^ stiffness only, a relatively simple relationship can be 

derived for the critical flutter parameter for panels vith relatively small 
D x stiffness. If Dq = oo, the constants in equation (25) become 
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TABLE II.- THE EFFECT OF TRANSVERSE SHEARING STIFFNESSES ON THE FLUTTER OF SIMPLY SUPPORTED 


PANELS FOR AN EXAMPLE SERIES (TWO-MOHE GAT.ERKTN SOLUTION) 





\:r 




>v cr (“) 

-A(r) 

Eq. (26): | =■ 1; (ly = 0.15; l,r,B = 1,2,1; 


Eq. ( 26 ): | = 1; Hy - 0.15; i,r,s - 1,2,1; 


^ » 1000 : C v = 0 ; N Y =■ 0 

Eq. (30): C y = 0 

D y 

e 5 - = 0.001; C v = 0.1, N y = 0 
% y x 


Eq. (29): Cy = 0; i,r = 1,2; 

*x - 1 


100 

1.0000 

1.0000 

0.8314 

1.0000 

300 

1.0000 

1.0000 

.7938 

1.0000 

1 000 

1.0000 

1.0000 

.6929 

1.0000 

3 000 

1.0000 

1.0000 

.4987 

1.0000 

10 000 

1.0000 

1.0000 

.2228 

1.0000 



c x = 1 

X 10 " 1 * 


100 

0.9943 

0.9950 

0.8267 

0.9950 

300 

■ 9845 

.9852 

.7820 

.9852 

1 000 

.9512 

• 9524 

.6652 

.9523 

3 000 

.8690 

.8696 

.4597 

.8695 

10 000 

.6664 

.6667 

.2078 

.6666 


c x = i x io "5 


100 
300 
1 000 
5 000 
10 000 

0.9456 

.8641 

.6638 

.3993 

.1668 

0.9524 

.8696 

.6667 

.4000 

.1667 

0.7866 

.6909 

.4891 

.2663 

.1024 

0.9519 
.8691 
.666 3 
.3998 
.1666 



Cx - 1 

_2 

X 10 


100 

0.6390 

0.6667 

0.5326 

0.6635 

300 

.3935 

.4000 

.3222 

.3980 

1 000 

.1678 

.1667 

.1356 

.1658 

3 000 

.06357 

.06250 

.05112 

.06218 

10 000 

.02003 

.01961 

. 01607 

.01951 



Cx = 1 

X 10' 1 


100 

0.1752 

0.1667 

0.1429 

0.1587 

300 

.07252 

.06250 

.05876 

.05942 

1 000 

. 02374 

.01961 

.01919 

.01863 

3 000 

.00812 

.006622 

.00656 

.006292 

10 000 

.00246 

.00200 

.00199 

.00190 


Ki(j) 


Kg(j) 



+ 


s 


2 af D xy\ 
t>2 Dx J 


l (28 


4.2 a 2 D xy 


J s 


b 2 D x 


+ 


0 .2 4 
2j s 





s^a^ ®y D xy 

.4 

b6 D x D x 

+ J 

— 



•^y x~y / \ 

Now, since ^ and are both large quantities (for small values of J 

of the same order of magnitude, only those terms in equations ( 28 ) that are 
underlined are retained. Therefore, A cr can be written (for any two modes 

and arbitrary stiffnesses) as 



(29) 


lllllllllllllllllllll 




jA(i 2 - r 2 ) 
8ir 


i 2 + r 2 - A 

! + c s 2 ^ 

1 + ° x b 2 D* 


where A(i) has become 


s - ( i+ 5 sfH 


2 a 2 Pxy ' 

h 2 W y y 


This expression is extremely easy to evaluate in comparison with equa- 
tion (26). Also, if A is substituted into equation ( 29 ) directly, the effect 
of in-plane loading is seen to be linear. Thus, an increase in compressive 
loading lowers the critical dynamic -pres sure parameter, and an increase in 
tensile loading raises it. This equation also shows that the influence of R^. 

decreases for the higher modes, and the largest effect is for i and r equal 
to 1 and 2. Finally, the influence of R x is greater for panels with smaller 


values of C x ^or larger values of Dq\ . A plot 


of 


A, 


\jr 

T 


with A for two 


representative values of R x and a series of values of C x is shown in fig- 
ure 7* The values used for plotting figure 7 are also given in table II for 
easy reference. 


The validity of the approximation used to obtain equation ( 29 ) may be seen 
if a comparison is made with the example given in figure 6. Thus, if N x = 0, 

equation ( 29 ) can be rewritten in the following form (for any planform ratio 
and any two modes): 


\r _ 1 

^cr^ 00 ) , AC X 

2~ 


(30) 


The dashed curves drawn in figure 6(a) are the flutter curves specified by 
equation (30). These values are also presented in table II. The approximation 
is extremely good in all regions and indicates that for this case, expres- 
sion (29) is accurate. It must be mentioned, however, that these conclusions 
are all based upon a two-mode Galerkin solution and it is known (ref. l) that 
the Galerkin two-mode solution is conservative when compared with the exact 
solution, particularly for large negative A values. 
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-100 -1000 - 10,000 


A 

Figure 'J.- The combined effect of Dq x and N x on the flutter of a simply 
supported orthotropic panel (two-mode Galerkin solution, Cy = o) for 
relatively small values of D x . 


CONCLUDING REMARKS 


Some important conclusions can be made on the effects of both si de-edge 
boundary conditions and transverse shearing stiffnesses on panel flutter phe- 
nomena. Since these two effects were investigated independently, separate con- 
clusions are presented. However, it must be emphasized that comparisons are 
made on the basis of two-term Galerkin solutions; exact solutions may not yield 
the same quantitative results, although the qualitative results are expected to 
remain the same . 
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Side-Edge Boundary Conditions 


The influence of side-edge boundary conditions can be very significant. 
Changing boundary conditions from simply supported to clamped causes a rela- 
tively small increase in the critical flutter parameter \» r ^ but changing 

these boundary conditions to free edges causes a significant decrease in A cr . 

This decrease is primarily caused (for the examples analyzed) by the fact that 
the panels were not forced to band between their side edges (their stiffest 
direction). Thus, if free boundary conditions exist along the side edges of a 
panel, flutter of that panel may occur at much lower dynamic pressures than 
expected for simply supported panels. 


Transverse Shearing Stiffnesses 

The two stiffnesses ^Dq x and Dq^ can cause extremely large variation in 

the critical flutter parameter some reductions being as large as two 

orders of magnitude. These conclusions are based upon flutter calculations for 
an example series of panels that closely approximate typical corrugation- 
stiffened panels currently in use, and for such panels (with relatively small 
bending stiffness in the flow direction), it was noted that if there is a 
finite transverse shearing stiffness in the flow direction ^Dq^, the effect 

of Dq^ is small. Therefore, a simple expression was derived for infinite 

Dq that showed the influence of Dq^_, in-plane load, and the mode numbers on 

the critical flutter parameter. This expression is extremely easy to use and 
agrees well with the more exact expression. Therefore, an est im ate of the 
flutter parameter may be easily calculated. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station,. Hampton, Va., September 17 , 19 & 5 . 
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APPENDIX 


GENERAL GALERKIN SOLUTION FOR THE FLUTTER OF ORTHOTROPIC PANELS 
IN SUPERSONIC FLOW USING STATIC STRIP THEORY 


The orthotropic panels under consideration here are shown in figure 1, 
with arbitrary boundary conditions on all four edges. The differential equa- 
tion for the lateral deflection w, both static strip theory and simple har- 
monic motion being assumed, is (eq. (l)) 


dNf + o/ D xy + \ S^w + c£w + fl2R x d 2 w + ^y <£w + _A_ dw 

\ D x / Sx 2 cly 2 D x Sy^ a^ dx 2 a 2 ay 2 s? ^ x 


£ 2 iA 


w = 0 


where 


(Al) 


N x a‘ 


R v = „ 

x ,2 


it D. 


N v a c 


A = 


2qa' 


-2 moi^a^ 
k = —i 


For the Galerkin solution, the deflection w was expanded in the series 



n p 


(A2) 


where each term ap n G n p(x) Fp n (y) must be chosen such that all the boundary 

conditions are satisfied. Substituting this function for w into equa- 
tion (Al) , multiplying by F^ m (y), and integrating across the panel yields (for 
each value of q) 
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y y bl( G )(m,n) G^(x) + — ( + y. \l( 2 )(m,n) G (x) + ^ P - n P^ X ^ I^^m^n) 

Z_. U qp ' pp b \d x y) qp np' d x b3 qp ' 


n p 


bR** 2 (0) 


V?« ^ _£2A» 


,2 ® 


Iqu ( m > n ) G nn( x ) + A- G nn( x ) I qb ( m > n ) 


"np'~' ’ o np'~' J -qp 
ba 


I — J qp A n ) G np( x ) 


. bA _(0). , ' , , 

+ — Iqp (m,n) G np (x) 

A 


^pn 


= 0 


(A3) 


for q equal to all values taken by p and where primes denote differentia- 
tion with respect to x and 



(m,n) 



■^qm(y ^ 


r pn(y) V 

dy 


Multiplying equation (A3) by G m q(x) and integrating along the panel yields 
(for each combination of q and m) 



p n <- 


J^U/P) I^(m,n) + ~ f=^ + 4 V^(q,p) I^(m,n) 

mn qp v ’ 1 b 2 \D X y mn qp v ’ ’ 


+ ^ Sy J^°\q,p) I^(m,n) + R n 2 J^ 2 ^(q,p) I^°^(m,n) 
D mn qp x ’ ' x mn qp v ’ ' 

X D 


2 ~ 2 (o), 


— 3f~a“ T ^ u >/ \ T ( 2 )/ \ . A 1 )/ \ T (°)/ \ 

+ *y ~ 2 ~ J mn (4,P) Iqp (^n) + AJ^ (q,p) Igp Kn) 


-2 4 t ( 0), v T (0)/ , 

k n J (q.p) I' (m.n) 
mn ' qp ' ’ ' 


a pn = 0 


(A4) 


for q equal to all values taken by p and for m equal to all values taken 
by n and where 
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mn 


(q>p) 



d r 

dx r 


G (x) 
np 4 


dx 


For any given boundary conditions on all four edges of the panel, func- 
tions are selected that satisfy all the boundary conditions. These functions 

are then used to calculate the integral expressions I^^(m,n) and J^^(q.,p) 

and the determinant of equation (Ak) yields a relationship to solve for the 
flutter parameter A. When generalized flutter parameters are introduced, 
equation (A4) becomes 


X 


*0 


<m 


T ( 0) / \ 

J (4) (q,p) + J (2) (q,p) 

mm (o), , mm qp 

I (m,m) 

qq ' 


J ( 0 ) (q,p) 

mm ' • gp 


mm 


I^(m,m) 

+ W^'(q,p) ® 

” I< 0 ) (m,m) 

qq ’ 


r(D/ 


XX 

P n 
n^m 


^n 


/ i \ I^°^(m,n) 

M, \ qp 

^mn ^q>P) 


P w (m,m) 

qq ’ 


(o), 


+ J ( 2 ) (q,p> K 2 A (m < n) - J (0) (q,p) A ( "’ n > + M (l) (q,p) ^ — 

mn qp mn qp mn (q^ 


I^°^(m,m) 
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= 0 (A5) 


for q equal to all values taken by p and for m equal to all values taken 
by n and where 
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If only one term is considered in the y-direction (as is often done), these 
equations become (for each value of m) 
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(A 6) 


for m equal to all values taken by n. 

As an example of the use of equation (A6 ) , the flutter parameter for a 
two-mode Galerkin solution can be written directly as a solution of the alge- 
braic equation 
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(Equation continued on next page) 
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These solutions are utilized in the text of this paper where applicable. 
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